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ARO  Final  Report 

Numerical  Modeling  of  Metamaterials 

Principal  Investigator:  N.  M.  Litchinitser 
Web  Designer:  V.  Z.  lakhnine 

Consultants:  V.  M.  Shalaev,  A.  V.  Kildishev 


1 .  Objective 

The  goal  of  this  project  was  to  build  a  website  containing  a  comprehensive  review  of 
available  software  packages  for  numerical  modeling,  design,  characterization,  and  applications 
of  metamaterials. 

2.  Summary 

We  have  developed  a  website  on  Numerical  Modeling  of  Metamaterials  (http://optical- 
waveguides-modeling.net/metamaterials/to  be  released  September  2010)  that  contains: 

•  Metamaterials  Tutorial  (as  illustrated  in  Fig.  1),  which  summarizes  basic  concepts  of  light 
propagation  in  optical  metamaterials,  magnetic  and  negative  index  metamaterials,  major 
challenges  and  approaches  to  numerical  modeling  of  metamaterials,  including  numerical 
studies  of  nonlinear  metamaterials,  multiscale  modeling,  materials'  parameters  retrieval 
procedures,  optimization  algorithms  and  a  summary  of  proposed  applications  of 
metamaterials; 
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Metamaterials  Tutorial:  Introduction 

Optical  Metamaterials 

Nobel  Laureate  Richard  Feynman  said,  "...The  imagination  of  nature  is  far, 
far  greater  than  the  imagination  of  man."  The  current  progress  in  the  field 
of  electromagnetic  metamaterials  ("meta"  stands  for  "beyond"  in  Greek) 
seems  to  challenge  this  statement. 

Despite  the  seemingly  unlimited  variety  of  physical  properties  of  the 
materials  existing  in  nature,  it  turns  out  that  only  a  very  limited  portion  of 
the  so-called  parameter  space  is  being  actually  utilized.  By  the  parameter 
space  we  understand  various  combinations  of  the  two  basic  parameters 
that  characterize  any  electromagnetic  medium  -  the  relative  dielectric 
permittivity  £  and  the  relative  magnetic  permeability  p,  or  their  product,  the 
index  of  refraction  n,  related  to  the  £  and  p  through  n  =  Vep 

To  start  with,  a  majority  of  optical  materials  such  as  water  or 
semiconductors  possess  £  that  is  greater  than  1.  Dielectric  permittivities  of 
metals  may  be  less  than  one  or  negative,  but  metals  are  lossy  in  the 
optical  range.  Moreover,  p  is  usually  close  to  one  at  optical  frequencies. 
Thus,  the  light-matter  interaction  is  dominated  by  the  interaction  with  the 
electric  field  component  of  the  electromagnetic  wave,  while  the  magnetic 
field  component  is  not  involved.  Quantitatively,  the  ratio  between  the 
magnetic  and  electric  coupling  strengths  is  of  the  order  of  10'2.  This 
significant  disproportion  stems  from  the  basic  properties  of  the  constituent 
components  of  any  naturally  existing  materials  -  the  atoms.  To  sum  up, 
optically  transparent  materials  are  almost  non-magnetic  and  possess  a 
positive  dielectric  permittivity  and  positive  refractive  index  that  is  typically 
greater  than  one.  But  is  there  something  unique  about  this  particular 
combination  of  £  and  p? 

Magnetic  or  even  negative 

In  1968,  Victor  Veselago  theoretically  investigated  a  class  of  materials  with 
both  £  and  p  being  negative  in  the  same  frequency  range  and  predicted 
that  such  materials  would  be  optically  transparent  and  would  give  rise  to  a 
number  of  extraordinary  phenomena  {Soviet  Physics  Uspekhi  10  (1968) 
509).  For  instance,  one  of  the  most  remarkable  properties  of  such 
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Figure  1.  Metamaterials  Tutorial  page 


•  A  broad  collection  of  the  software  available  for  numerical  simulations  of  metamaterials 
structures,  supplied  with  a  short  summary  of  its  capabilities  and  potential  applications,  a  list 
of  references  to  research  papers  that  utilize  a  particular  software  package,  and  a  link  to  the 
software  provider's  page. 

•  We  provide  the  visitors  of  our  website  with  an  online  file-sharing  facility  to  be  used  to 
exchange  simulation  codes,  documentation,  and  other  relevant  information.  This  facility  is  a 
first  step  toward  building  an  online  research  community  with  its  members  sharing  news, 
scientific  expertise,  and  knowledge  in  the  rapidly  developing  field  of  optical  metamaterials 
modeling. 

•  The  website  contains  a  number  of  simulations  demos  as  well  as  a  link  to  a  number  of 
simulation  tools  available  on  nanoHUB.org 

•  Near  term  future  plans  include  creating  a  discussion  group  and  an  online  forum. 

3.  Website  Structure  Overview 

The  website  contains  the  following  pages:  Home,  Metamaterials  Tutorial,  Software, 
Bibliography,  Site  Map,  Community,  and  Search.  The  front  page  contains  a  “Welcome” 
message  and  a  brief  introduction  to  the  website  structure.  It  also  includes  a  “News”  section  with 
announcements  of  recent  breakthroughs  in  the  field  of  metamaterials  and  a  list  of  upcoming 
research  meetings,  workshops,  and  conferences  related  to  metamaterials  and  plasmonics. 
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Optical  Metamaterials: 

1 

Numerical  Modeling 

Home  |  Metamaterials  Tutorial  4-  Software 

4>  Bibliography  Site  Map  |  Community  4>  |  Search  | 

VECTOR  FIELDS:  CONCERTO 

(-fURi  mm  Concerto  is  an  advanced  analysis  package  for  RF 
and  Microwave  design.  It  includes  a  3D  Geometric 
Modeler  that  allows  easy  data  import  from  CAD 
systems,  as  well  as  the  facility  to  create  and 
WSW^M  parameterize  new  models. 

) 

General  features 
Typical  applications  include: 


•  Metamaterials 

•  Cavities 

•  Antenna  design 

•  Microwave  component  design, 
including: 

■  Filters 

•  Diplexers 

•  Couplers 

•  Junctions 

•  Phase  shifters 

•  Polarizers 

•  Microstrip  component  design 

•  Microstrip  to  waveguide  transitions 

Concerto  includes  the  following  modules: 


B.  Hou  et  al.,  Tuning  of  photonic 
bandgaps  by  a  field-induced 
structural  change  of  fractal 
metamaterials,  Opt.  Express  13. 
9149-9154  12005). 


>  Concerto  AS:  Introducing  Low  Cost  Solution  for  Coaxial  Connector 
Design. 

The  FDTD  method  is  used  for  the  specific  task  of  designing  coax 
connectors  using  a  fast  and  efficient  2D  axi-symmetric  analysis. 

>  Concerto  ES:  Entry  Level  3D  FDTD  analysis. 

The  Concerto  ES  module  includes  the  highly  efficient  FDTD  method,  with 
options  to  add  other  analysis  modules.  The  entry  price  is  kept  very 
competitive  by  restricting  the  size  of  model  that  can  be  solved,  whilst 
retaining  its  ability  to  analyze  a  wide  range  of  devices. 

>  Concerto  2D:  Axi-symmetric  Analysis  using  2D  FDTD. 
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Figure  2.  An  example  of  software  description  page  (providing  a  general  description  of  the  software,  typical  applications 

and  basic  modules). 


The  Software  page  is  a  key  part  of  the  website.  It  contains  an  extensive  collection  of  software 
packages  for  various  types  of  metamaterials  simulations,  including  linear  and  nonlinear  wave 
propagation,  frequency  and  time  domain  methods,  including  finite  difference  time  domain 
method,  finite  elements  method,  method  of  moments,  and  eigenmode  solvers,  a  description  of 
typical  applications  of  a  particular  software  package,  hardware  requirements,  and  some  useful 
references,  as  shown  in  Figs.  2  and  3. 

For  users  who  are  new  to  the  field  of  metamaterials,  we  offer  a  short  Metamaterials  Tutorial, 
which  summarizes  basic  concepts  of  light  propagation  in  metamaterials,  metamaterials  design, 
optimization,  and  characterization. 

We  compiled  a  Bibliography  that  contains  many  useful  references  about  metamaterials, 
plasmonics,  nonlinear  optics,  and  other  general  references  on  nanophotonics  and 
nanotechnology. 

The  website  is  also  supplied  with  a  Site  Map  for  easier  navigation  through  the  site,  a 
“community”  section,  providing  an  online  file-sharing  facility  to  be  used  to  exchange  simulation 
codes,  documentation,  and  other  relevant  information,  as  well  as  a  search  engine,  enabling 
visitors  to  explore  the  site  with  the  keywords  of  their  choice. 

In  addition  to  the  detailed  description  of  the  software  available  on  the  web,  the  website 
contains  a  summary  of  modeling  tools  developed  largely  by  Purdue  University  researchers  and 
students  that  is  available  on  Nanohub,  as  shown  in  Fig.  4.  We  summarize  the  basic  features  of 
the  tools  related  to  numerical  modeling  of  metamaterials  and  plasmonic  nanostructures,  and 
provide  a  link  to  the  tool  itself.  One  of  the  scientific  consultants  involved  in  this  project,  Dr.  A.  V. 
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=  Ideal  for  antenna  coupling,  antenna  installed  performance 

=  Radar  signature  calculations  (including  mono-static  and  bi-static 
RCS) 

•  Customized  post-processor  for  result  display  and  further 
processing 

Finite  Element  Method 

=  Eigenvalue  analysis  of  resonant  structures  analysis  (single 
frequencies) 

--  Uses  finite  element  method  for  efficient  modeling,  with  automatic 
mesh  generation 

=  Ideal  for  RF  and  other  cavities 

»  Customized  post-processor  for  result  display  and  further 
processing 


1.  X.  Huang  et  al.,  Fractal  plasmonic  metamaterials  for  subwavelength 
imaging,  Opt.  Express  18.  10377-10387  (2010). 

2.  H.  Xu  et  al.,  Effective-medium  models  and  experiments  for  extraordinary 
transmission  in  metamaterial-loaded  waveguides,  Appl.  Phys.  Lett.  92, 
041122  (2008). 

3.  3.  Hao  et  al.,  Optical  metamaterial  for  polarization  control,  Phys.  Rev.  A 
80,  023807  (2009). 

4.  H.  Li  et  al.,  All-dimensional  subwavelength  cavities  made  with 
metamaterials,  Appl.  Phys.  Lett.  89,  104101  (2006). 

5.  L.  Zhou  et  al.,  Electromagnetic-Wave  Tunneling  Through  Negative- 
Permittivity  Media  with  High  Magnetic  Fields,  Phys.  Rev.  Lett.  94,  243905 
(2005). 
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Figure  3.  An  example  of  software  description  page  (providing  a  description  of  numerical  methods  implemented  in  this 
particular  software  package,  a  summary  of  different  features,  and  some  useful  references). 
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Nanohub  Tools 

PhotonicsSHA-2D:  Modeling  of  Single-Period 
Multilayer  Optical  Gratings  and  Metamaterials 

PhotonicsSHA-2D: 

Modeling  of  Single-Period 
Multilayer  Optical  Gratings 
and  Metamaterials 
employs  the  2- 
dimensional  spatial 
harmonic  analysis  (2D 
SHA)  method,  also  known 
as  Fourier  modal  method 
(FMM)  and  rigorous 

coupled-wave  analysis  (RCWA),  to  model  the  electromagnetic  behaviors  of 
single-period  multilayer  gratings.  The  incident  wave  is  a  plane  wave  with 
arbitrary  incident  angles  (with  either  TE  or  TM  polarization),  and  the  output 
results  are  the  complex  transmission  and  reflection  coefficients  for  the  zero 
diffraction  order.  The  database  of  optical  elemental  materials  stored  in 
PhotonicsDB  is  integrated  into  this  tool.  An  online  video  tutorial  is  also 
available.  To  learn  more  and  try  the  tool  out,  follow  this  link. 

Hyperlens  Design  Solver 

Recent  research 
has  been  done 
in  regard  to 
optical  imaging 
using 

metamaterials. 

One  such  project 
is  the  hyperlens, 
which  aims  to 
overcome  the 

classical  diffraction  limit  and  project  a  magnified  image  into  the  far  field.  The 
potential  applications  for  this  device  range  from  nanolithography  to 
bioimaging.  The  Hyperlens  Design  Solver  tool  is  intended  to  be  used  in 
conjunction  with  the  Hyperiens  Layer  Designer  tool  to  aid  in  the  design  and 
simulation  of  a  hyperlens.  The  Hyperlens  Design  Solver  tool  allows  users  to 
upload  designs  created  with  the  Hyperlens  Layer  Designer  tool,  make  a 
new  design,  or  select  from  several  pre-existing  designs.  The  tool  then 
simulates  the  performance  of  the  design  and  outputs  several  plots  of  the 
resulting  field  intensities.  By  using  these  two  tools,  users  can  experiment 
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Figure  4.  A  summary  of  Nanohub  tools  with  links. 

Kildishev,  has  developed  a  number  of  tools  available  on  Nanohub,  and  therefore,  our  website 
users  will  be  able  to  obtain  a  detailed  advice  and  guidance  regarding  those  tools,  if  needed. 

Finally,  several  online  demos  have  been  developed  enabling  our  users  and  visitors  to  watch 
demo-simulations  of  such  fascinating  metamaterials-based  devices  as  a  cloak  or  a  hyperlens, 


Figure  5.  Examples  of  online  demos  for  a  cloak  (left)  and  a  hyperlens  (right) 


as  shown  in  Fig.  5. 


The  metamaterials  modeling  website  will  be  continuously  maintained,  including  regular  content 
updates  (e.g.,  adding  new  tools  as  they  become  available),  small-to-medium  page  redesigns 
(e.g.,  incorporation  of  a  videostream,  a  slide  show,  or  another  media  type),  and  creation  of  new 
pages,  as  needed,  with  subsequent  modification  of  the  site  navigation.  Usage  statistics  of  the 
sites  will  be  monitored  and  analyzed,  as  well  as  sites’  ranking  by  major  search  engines. 
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